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a Have a mesh with its boundary elements. Pick one. 
a Heat as blackbody radiation radiates from this face 
a Heats up other faces the radiation impinges 


A view factor quantifies radiation transfer between areas: 

radiation leaving Ai and impinging upon A2 


F( Ar,A 2 ) 


all radiation leaving Ai 
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Binary Trees, 

Quadtrees, 

and Octrees 



a Can compute view factors by Monte Carlo ray casting 
a No info in ray on where it lands; must check all faces in set 
X Ok for small meshes, but gets glacially slow quickly 
a Binary trees used in computer science for efficient searches 
a Because meshes are 2D and 3D, quadtrees and octrees are 
the weapon of choice to speed up the searches 
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“Glacially Slow”? 


Orion tile cavity problems 
take too long: 



“The Monster” 
26 , 000 + faces 

Bonus: add capability to solve 
Eric’s fiber problem. 
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A Tree Grows 

in Data 




Suppose a data set that has some sort of spatial association. 


Bin Size b: 
99 

32 

16 

8 

4 



Binary trees are formed by dividing data set in half repeatedly. 
Do that in each dimension until each box has b items. 

Data associated with space: Boundary faces T 
Placement criterion: T ^ B if \T n B\ d _ 1 > 0 
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The Hyperplane Separation Theorem 


We pull out the artillery: convex analysis. We will power our 
intersection algorithm with this theorem: 


Theorem (Hyperplane Separation Theorem) 


Let A c M. d be closed and /C c W d be compact with both convex. 
Then A n /C = 0 if and only if there is a separating hyperplane 
V = jx e : x • p = a} for some a € M and p e M d \{0} such that 

9 p • a > a for all p e V, a e A and p • k < a for all p e V, k e K. 

or 

Q pa < a for all p e V, a e A and p • k > a for all p e V, k € K, 


Theorem (Plain Language Version) 


Two nice sets A and /C don’t intersect if we can divide space 
into A’s half and JC 1 s half with a line (2D) or plane (3D). 
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The Separating Axis Theorem 


Two vectors span a plane in 3D, so testing for hyperplanes 
directly is expensive. Cheaper to test for separating axes: 


Definition (Separating Axis) 


Let V c M. d be a separating hyperplane. £ c is a separating 
axis if £ 1 V. 

Because dimV = d- 1, dim£ = 1. This leads to an obvious result: 


Theorem 


Let i c K d be closed and /C c be compact with both convex. 
Then there exists a separating hyperplane for A and K. if and 
only if there exists a separating axis between them. 

The key: when A and K are orthogonally projected onto 
A fj K, = 0 if and only if the projection intervals do not overlap. 
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Collecting the Candidates 




There are six ways for intersection to occur: 


» face-face 
o face-edge 
<» edge-edge 


Same 
in 2D 


y face-node 
edge-node 
y node-node 

Can perturb non-intersecting faces and still maintain 
separation. Gives these candidate separating axes: 

Q Face normals from /C 

Q Face normals from K' / 


K 


Q Cross product of edge from /C 
with one from K.' (3D only) 


K' 


In 2D, these are sufficient. In 3D, they cover face-face and 
face-edge cases. Edge-edge cases can look like this — 
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SAT for Symmetric Objects 


9 Let c*; € fC and c k ' e JC' be the centroids of K. and K! . 

Q Let n(x) denote the projection of x e W d onto Define 
r k = max |n(nj-c fc )| r k > = max |n(n--c*/)| 

p - |n(c fc - c fc OI 

O If r k + r k' < Pi then £ is a separating axis. 
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How It Works 


The SAT for symmetric objects works so efficiently because 
Q K. and 1C' are symmetric about centroids C& and cy. 

=> Projection intervals symmetric about n(cj,) and II(c' fc ). 
=> It is sufficient to project the “radii” of 1C and /C'. 

O By linearity of projection operator, we can project the 
centroid-to-centroid segment. 

For non-symmetric objects, it is a bit more complicated: 

Q Must project every node of K. and K! . 

Q Must project every node of the convex hull of /C and 1C' . 

Because 1C and 1C 1 are convex, extreme points of convex hull 
come from nodes of 1C and 1C' . Can reuse values from Step 1. 
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SAT for Non-symmetric Objects 
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Algorithm 


Q Let {n*}^ 1 and {n^]T 1 be nodes 
of /C and K.' . Put 


i JV'-l 


Pi 


n, 


p'i - K A 


Compute lengths of proj. intervals: 


n 


max pi 
0<i<N—l 


mm pi 
0<i<N—l 


ruf = max p [ - min p- 
0<i<N-l 0<i<N-l 


max < 


max pi , max 
0<j<iV-l 0 <i<N- 


/i) 


mm 


in{ min Pi , min p'A 
( 0<i<N-l 0<i<N-l J 


0 £ is separating axis if r & + ry < p ■ 
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Candidates for separating axes are face normals of box and 
normal to the face/edge w = ni - no: 

£i = (i,0) T £ 2 = (0,1) T £ 3 = (-™ (1) ^ (0) ) t 


Both B and T are symmetric: can use SAT for symmetric 
objects. Project w, box radius d = x max - x m j n , and 
centroid-to-centroid vector m = no + ni - x max - x m ; n onto £: 


Axis 

€ 

rk 

rk' 

P 

1 

(1,0) T 

\d (0) \ 

l^°)| 


2 

(0,1) T 


]«;«] 

|m«| 

3 

1 

o 

i d ( o ) w ( i )i + i d ( i ) w (°)i 

0 



Because > 0, we can lose some absolute values in column 3. 
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In 3D, we do not expect T to be symmetric. But B is! Can 
compute the max and min of projection nodes Vj of B directly: 

M = max(vj • £) - £ ' x + m= min(vj • £) = ^ • x_ 


| XmL if > o = | if e (i) > o 

USn if C (i) < 0 X% U^L if£ (i) <0 


First set of candidate are B’s face normals = (1,0, 0) 1 ", 
= (0, 1,0) T , and £ 3 = (0,0, 1) T and normal v from T\ 


Axis 

£ 

rk 


p 

i 

(1,0, 0) T 

d(°> 

(0) . (0) 
max n\ - mm n ^ 

( (0) (0)1 • f (0) • (0)1 

max \Xmax , max n\ mm \ x ^{ n , mm n) ) 

2 

(0,1, 0) T 

dW 

(i) • (i) 

max - mm 

r (i) (l)i • r (l) • (i)i 

max |#max? max \ - mm uV n , mm n\ \ 

3 

(0,0, 1) T 

d< 2 > 

(2) . (2) 
max - mm 

f (2) (2) l . r (2) . (2) l 

max jaimax, max j - mm j x ^„{ n , mm n\ j 

4 

V 

M-m 

max v ■ rii - min v ■ n ? 

max {Af, max v ■ n^} - min {m, mini/ • n^} 
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SAT for Insertion in 3D (continued) 


Second set of axes given by cross products of box edges with 


face edges f \ = n, + i - n,: 


£5,1 = (1)0, 0) T x fj 

a j,i - n j ' £5,1 

£ 6ii = (o,i,o) T xfc 

Pj,i = n j ' ^ 6 ,i 

| 7jJ = (0,0,l) T xf ? ; 

7 J,i ~ n j ' %7,i 


Can compute everything without forming f \ or ever computing 
cross product. Also have = ctj+i,i, similarly for and 7^. 


Axis 

€ 

rk 

rk 1 

p 

5,2 

h ^ 

cu 

1 

O 

A/5 — 77 l 5 ,t 

max aj i - min a 7 - i 
j*i j*i 

max {Af5 )7 , maxajd - min {7715,4, min ocj t i] 
j*i ’ ’ j*i 

6, i 

O 

1 

0 

A/(6 ,i — 2716 , z 

max - min fijj 

j*i ’ 3 ±i 

max { Mqj , max ( 3 j * } - min { rri 6 ,i , min f 3 jj } 
j±i j±i 


1 

0 

0 

4 

- mr,i 

max 7 jj, - min 7 jj 

jti jti 

max { M7 * , max 7 ji } - min { 7717,4 , min 7 ji } 
j*i j*i 


Mjj and rrij ^ are dehned/computed like M and m of last slide. 
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Trees Are Made for Climbing 


a Have a mesh. With a tree built from it. 
a Shoot a ray from one of the faces, 
a Now want to identify which boxes the ray visits. 
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Intersection with Child Boxes 


SAT gives no info about intersection point, so poor for 
recursion. Traversal algorithm will look like “first attempt.” 

Suppose we have ray t(f ) = p + tr. p is inside mother box, so we 
need to move it out: p p - nr At, where 

„(0 1 


At = - min 

5 1 <i<d 


' (<) 
*^max 


XI 



I |r(*)| j 
Compute times: 


Entry : 


(i) ( i ) 

x entry ”0 


Exit 
Midpoint : t. 


7"(*) 
+(*) ‘^exit 

h “ M 


(i) (i) 

xLh - V o 


(0 A 
m 


r (i) 

lew , ,(0' 
2 L 0 1 
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2D Ray Traversal 



Child Box 

Condition for Entry 

Entry Times 

Exit Times 

vo 

max 

KW) 

< min < 

.(0) ,(1)1 
*'777 5 *'777 [ 

,(0) 

r 0 

t (1) 

l 0 

,(0) 

*-771 

t (1) 

*'771 

Vl 

max 

{4°DT} 

< min < 

.(0) ,(i)l 

5 [ 

AO) 

z o 

t w 

*'771 

AO) 

*'771 

t (1) 

V2 

max 

/ t <°> t (1) ) 

1 *'771 5 ‘'O I 

< min < 

.(o) f W\ 

‘'l 5 *'771 j 

A0) 

L m 

l 0 

f (0) 

6 1 

t (1) 

*'771 

V3 

max 

h«» t (1) ) 

| *'771 5 *'771 

< min < 

4VT} 

A0) 

'-"m 

f (D 

f (0) 

t (1) 

6 1 


...but this table will only be valid if r h) > 0 for each i. 
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2D Ray Traversal 


Time spent in: 



Child Box 

Condition for Entry 

Entry Times 

Exit Times 

ro 

max 

KW) 

< min < 

.(0) ,(1)1 

,(0) 

l 0 

t (1) 

l 0 

Ao) 

^m 

t (1) 

L m 

Vl 

max 

{4°DT} 

< min < 

.(0) ,(1)1 
^ m 5 [ 

AO) 

z o 

t w 

AO) 

t/ m 

t (1) 

L 1 

V2 

max 

/ t <°> t (1) ) 

< min < 

.(0) fW\ 

? L m j 

Ao) 

L m 

t (1) 

l 0 

Ao) 

L i 

t (1) 

v m 

V3 

max 

h(°> t (1) ) 

j b m i 

< min < 

4VT} 

Ao) 

l m 

f (D 

AO) 

L i 

t (1) 

6 1 


...but this table will only be valid if r 0) > 0 for each i. 


x m 

.mid 
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Rays with Negative Components 


Consider the reversal of the previous example: 



A simple relabeling of boxes will allow reuse of previous table! 
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The Relabeling Scheme 


The general mapping for box relabeling is £ i ® a, where 

A [ 2<r(0) + cr(l) if d = 2 

| 4<r(0) + 2cr(l) + c(2) if d = 3 

... f 0 if j-b) > 0 
= { 1 if rb) < 0 


Written as binary string, a encodes the signs of the components 
of r (1 negative, 0 nonnegative). 

In 3D, the box labeling order 
compatible with the relabeling 
scheme is this goofy thing: 


With this, we can make a table 
for 3D very similar to 2D one. 
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Performance 


One Last Thing: Infinite Arithmetic Module 


a Computation of entry/exit times totally fails if r ^ = 0. 
a If | A') | < e « 1, then the numerics are bad too. 
a Small (positive) direction value is inducer of exit only if 




> (x^ 

Vma? 




2e 2 


Solution: when | r ^ [ < r, set 

.(i) 

T v 7 — 


*(0 


+ oo 


a Says that xO n < t^\t) < x fiL for all t. 
a Computation of as before is undefined. Instead, define 


go 


+oo 


if ph) < x 
if u (i) > 


(0 

mid 

(0 

mid 
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768 faces 

i 

192 faces 


3D Concentric Spheres (5k Rays) 

768 faces, Trunk: 47.69s 


Bin Size 

Run Time 

Speed Up 

Rays Diff. 

Corrected 

30 

9.72s 

4.9x 1 



20 

7.99s 

6. Ox 

L 2452 

0.063854% 

10 

5.97s 

8. Ox 

5 

4.44s 

10. 7x J 




3D Concentric Spheres, 2-Plane Symmetry (5k Rays) 

192 faces, Trunk: 20.52s 


Bin Size 

Run Time 

Speed Up Rays Diff. 

Corrected 

30 

2.59s 

7.9x \ 


20 

2.15s 

l 2951 

0.255420% 

10 

1.86s 

11. Ox | 


5 

1.51s 

13. 6x J 
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764 faces 



1528 faces 


3D Cylinder, 1-Plane Symmetry (10k Rays) 

764 faces, Trunk: 77.22s 


Bin Size 

Run Time 

Speed Up Rays Diff. 

Corrected 

30 

17.16s 

4.5x \ 


20 

9.03s 

* 

0.000104% 

10 

7.37s 

10. 5x J 


8 

5.90s 

13. lx J 



3D Cylinder (10k Rays) 

1528 faces, Trunk: 307.05s 


Bin Size 

Run Time 

Speed Up Rays Diff. 

Corrected 

30 

41.41s 

7.4x 1 


20 

21.83s 

14. lx 




^ 16 

0.000104% 

10 

17.95s 

17. lx | 


8 

14.07s 

21. 8x J 
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6414 faces 


n 

216 faces 



3D Cylinder, Fine Mesh (10k Rays) 



6414 faces, Trunk: 9892.82s 


Bin Size 

Run Time 

Speed Up 

Rays Diff. 

Corrected 

30 

203.97s 

48. 5x 



20 

10 

121.03s 

84.40s 

81. 5x 
117. 2x 

34 

0.000053% 

8 

75.02s 

131. 9x 




2D Cavity (5k Rays) 



216 faces, Trunk: 

3.01s 


Bin Size 

Run Time 

Speed Up 

Rays Diff. 

Corrected 

30 

1.48s 

2. 0x 



20 

1.39s 

2.2x 

° 


10 

1.27s 

2.4x 


5 

1.02s 

3. Ox 
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‘The Monster’ 



The primary motivation for doing all this work was this 

26,232 Faces 
1,000,000 Rays 


Force” 

96 hours 
12 Restarts 


Tree 

31 minutes 
0 Restarts 
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Summary and Conclusions 




Tree-based search: 

» Has three speeds: Fast, Blazing, and Ludicrous 
a Numerically robust thanks to Separating Axis Test 
a Pretty simple to implement 
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